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Abstract — This paper addresses the problem of controlling 
a variable-speed wind turbine with a Doubly Fed Induction 
Generator (DFIG), modeled as an electromechanically-coupled 
nonlinear system with rotor voltages and blade pitch angle as 
its inputs, active and reactive powers as its outputs, and most of 
the aerodynamic and mechanical parameters as its uncertainties. 
Using a blend of linear and nonlinear control strategies (including 
feedback linearization, pole placement, uncertainty estimation, 
and gradient-based potential function minimization) as well as 
time-scale separation in the dynamics, we develop a controller 
that is capable of maximizing the active power in the Maximum 
Power Tracking (MPT) mode, regulating the active power in the 
Power Regulation (PR) mode, seamlessly switching between the 
two modes, and simultaneously adjusting the reactive power 
to achieve a desired power factor. The controller consists of 
four cascaded components, uses realistic feedback signals, and 
operates without knowledge of the Cp-surface, air density, friction 
coefficient, and wind speed. Finally, we show the effectiveness of 
the controller via simulation with a realistic wind profile. 

Index Terms — Wind energy, wind turbine, active power, reac- 
tive power, maximum power tracking, power regulation, nonlin- 
ear control. 



I. Introduction 

A high-performance controller is essential to the success 
of integrating large-scale wind energy into future power 
systems. Extensive investigations and recent lessons learned 
have confirmed that the variable and intermittent nature of 
wind indeed poses serious threats to both the reliability of 
power systems and the economic viability of wind energy [JJ. 
To minimize these threats, a wind turbine controller should 
not only maximize the amount of active power captured 
in a so-called Maximum Power Tracking (MPT) mode in 
normal situations, it should also allow the power captured be 
continuously regulated at a desired level in a so-called Power 
Regulation (PR) mode (other than clipping the power based 
on received instructions) when there is a system contingency. 
In addition, the controller should enable seamless switching 
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between the MPT and PR modes, as well as maintain a desired 
power factor by also controlling the reactive power output. 

One of the challenges facing the development of such a 
high-performance controller is the fact that the aerodynamic 
and mechanical parameters of a wind turbine are inherently 
uncertain, due to modeling and measurement errors, unknown 
optimal operating points, and other, possibly time-varying, 
ambient factors. For example, the Cp-surface of a wind 
turbine, which characterizes the amount of mechanical energy 
converted from the wind, is typically assumed to be known — 
or, at least, its optimal points are assumed to be known — in 
many existing controller designs. Unfortunately, such a surface 
is an empirical, statistical approximation, obtained based on 
up to three months of continuous experiment f2\. Thus, the 
Cp-surface may not be precisely known for control purposes. 
Other factors, such as changes in air density due to weather, 
variations in friction under different operating conditions, and 
measurement errors due to anemometer location, also con- 
tribute to the uncertainties. Indeed, a report from the National 
Renewable Energy Laboratory (NREL) |3|, which describes 
the result of its long-term test on different controllers operating 
on real, MW-level wind turbines, shows that the impact of 
these uncertainties on controller performance is significant and 
should be accounted for in controller design. 

Another major challenge facing the development of such 
a high-performance controller is the fact that the mechanical 
and electrical parts of a modern wind turbine, which uses a 
Doubly Fed Induction Generator (DFIG), are tightly coupled. 
As will be detailed below, most studies have adopted a 
standard approach in the analysis and control of synchronous 
electric machines, in which the active and reactive powers 
are considered decoupled. With this approach, the active and 
reactive powers are adjusted via control of the mechanical 
and electrical parts, respectively, independent of each other. 
However, although a DFIG has some features of a synchronous 
machine, it is by nature an induction machine that exhibits 
strong electromechanical coupling among its rotor excitation 
current, rotor angular velocity, and electromagnetic torque. 
Hence, for performance reasons, the mechanical and electrical 
parts of a wind turbine with a DFIG should be considered 
synergistic ally in controller design. 

The current literature offers a large collection of wind 
turbine controllers, including Il4l- ll24l . However, as was al- 
luded to above, most of the existing publications considered 
the mechanical and electrical parts separately (e.g., ll4l- lfT4ll 
considered only the former, while flSl - llSSl considered only 
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the latter), and for a few of those (e.g., 112411 ') that considered 
both parts, its controller can only operate in the MPT mode, 
maximizing wind energy conversion, as opposed to both the 
MPT and PR modes. Moreover, although the existing work 
has provided valuable understanding in the control of wind 
turbines, only a few publications have addressed the issue of 
uncertainties. For example, |l4l-(|6l proposed adaptive frame- 
works for controlling the mechanical part of wind turbines, 
so that the power captured is maximized, despite not knowing 
the Cp-surface. 

In our recent work f25\, we developed a nonlinear controller 
that simultaneously enables control of the active power in 
both the MPT and PR modes, seamless switching between 
the two, and control of the reactive power so that a desirable 
power factor is maintained. These objectives were achieved by 
adjusting the rotor voltages of the electrical part and the blade 
pitch angle of the mechanical part, where the coupling between 
the two parts were taken into account in the controller design. 
Like most of the existing work, however, the controller in 
[jSyi assumed that the aerodynamic and mechanical parameters 
were known. 

In this paper, we develop a controller that achieves such 
objectives and, at the same time, addresses the two aforemen- 
tioned challenges, on uncertainties in the aerodynamic and me- 
chanical parameters, and coupling between the mechanical and 
electrical parts. For the former, we show that the parametric 
uncertainties can be lumped into a scalar term, estimated via 
an uncertainty estimator in an inner loop, and circumvented in 
an outer, gradient-based minimization loop. For the latter, we 
show that the electromechanical coupling can be eliminated via 
feedback linearization on the electrical dynamics, following 
ideas from 125|. Finally, we demonstrate the effectiveness of 
the controller developed through simulation with a realistic 
wind profile from a wind farm in Oklahoma. 

The outhne of this paper is as follows: Section|ll]models the 
wind turbine and formulates the problem. Sectionllllldescribes 
the proposed controller Section |IV] presents the simulation 
results. Finally, Section |V] concludes the paper. The proof of 
the main theorem is included in the Appendix. 



II. Modeling and Problem Formulation 

Consider a variable-speed wind turbine with a Doubly Fed 
Induction Generator (DFIG). The wind turbine consists of an 
electrical part and a mechanical part, the dynamics of which 
may be modeled as follows: 

The dynamics of the electrical part in the dq frame are 
described by a fourth-order state space model 1261 . ll27l 
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where ip^s, Vqs^ Vdr, 'fiqr G ^ ^"6 State variables representing 
the stator and rotor fluxes, Vdr, £ R are control variables 
representing the rotor voltages, Vds, Vqs G M are the constant 
stator voltages (that are not simultaneously zero), > is 
the constant angular velocity of the synchronously rotating 
reference frame, ujr > is the rotor angular velocity, Rs, 
Rr are the stator and rotor resistances, Ls, Lr, Lm are the 
stator, rotor, and mutual inductances satisfying Ls > L„i and 
Lr > L„i, (7 = 1— ^ is the leakage coefficient, and A, B 
are constant matrices. In addition, the fluxes can be written as 
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M are the stator and rotor currents, 

ipqr]'^ and i = [ids iqs idr iqrT 

introduced just for convenience. Furthermore, the active and 
reactive stator and rotor powers are given by Ii28l 

Ps ^ds^ds ^qs^qs^ Qs — '^qs'^ds ~^ '^ds'^qsj (3) 

L*r — ^dr^dr '^qr^qr, Qr — ^qr^dr '^dr^qr, (4) 

and the total active and reactive powers of the turbine are 

P = Ps+Rr, Q = Qs+Qr, (5) 

where positive (negative) values of R and Q mean that the 
turbine injects power into (draws power from) the grid. 

The dynamics of the mechanical part are described by a 
first-order state space model 1261 



JuJr — Tm — Pe — 



(6) 



where the rotor angular velocity ujr is another state variable, 
J is the moment of inertia, C/ is the friction coefficient, Tm 
is the mechanical torque, and Te is the electromagnetic torque 
given by [28] 



Vqs'^ds Vds'^qs-; 



(7) 



where positive (negative) value of Tg means that the turbine 
acts as a generator (motor). The mechanical power captured 
by the wind turbine is ||29l 



Rm = Tm^r = ]^pACp{\,P)V^, 



(8) 



where p is the air density, A — ttR^ is the area swept by the 
rotor blades of radius R, Vw is the wind speed, and Cp(A, (3), 
commonly referred to as the Cp-surface, is the performance 
coefficient of the wind turbine, whose value is a function of 
the tip speed ratio A £ (0, oo), defined as 

A=^, (9) 
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and the blade pitch angle /? G [f3min, f^max], which is another 
control variable. 

In order for results of this paper to be applicable to a broad 
class of wind turbines, no specific expression of Cp(A, /3) will 
be assumed. Instead, Cp{X, (3) will only be assumed to satisfy 
the following mild conditions for the purpose of analysis: 
(Al) Function Cp{X, /3) is continuously differentiable in both 

A and 13 over A e (0, oo) and /3 € [/3min, /3max]- 
(A2) There exists c £ (0, oo) such that for all A S (0, oo) 
and (3 £ [Amn,/3max], wc havc Cp(A,/3) < cA. This 
condition is mild because it is equivalent to saying that 
the mechanical torque is bounded from above, since 



according to ([8]l and (|9|l. 



(A3) For each fixed /3 e [/3min, /3max], there exists Ai £ 
(0, (X)) such that for all A € (0, Ai), we have Cp(A, /3) > 
0. This condition is also mild because turbines are 
designed to capture wind power over a wide range of 
A, including times when A is small. 

(A4) There exist c G (— oo, 0) and c E (0, oo) such that for 
all A e (0, oo) and 13 e [/3mm, 13 max], we have c < 

fCAMl] < r 
d\\ A I - 

As it follows from the above, the wind turbine is modeled 
as a fifth-order, electromechanically-coupled, nonlinear system 
with state equations (HJ and (|6]l, output equations (|3]l-(|5]l, state 
variables ipds, 'figs, 'fidr, 'Pgr, and ojr, control variables Vdr, Vqr, 
and /3, output variables P and Q, and exogenous "disturbance" 
Vw ■ A block diagram of this system is shown on the right-hand 
side of Figure [T] in which the electromechanical coupling can 
be seen. 

Given the above model, the problem addressed in this paper 
is: design a feedback controller, so that the active and reactive 
powers P and Q closely track some desired, possibly time- 
varying references Pd and Qd, assumed to be provided by a 
wind farm operator When Pd is larger than what the wind 
turbine is capable of generating, it means that the operator 
wants the turbine to operate in the Maximum Power Tracking 
(MPT) mode; otherwise, the Power Regulation (PR) mode is 
sought. By also providing Qd, the operator indirectly specifies 
a desired power factor PF^ = —^^i==, around which the 



actual power factor PF 



should be regulated. The 



controller may use i, ojr, P, and Q, which are all measurable, 
as feedback. The fluxes ip may also be viewed as feedback, 
since they are bijectively related to i through (|2]i. Moreover, 
the controller may use values of all the electrical parameters 
(i.e.. Us, Rs, Rr, Ls, Lr, L,n, Vds, and Vqs) and turbine- 
geometry-dependent parameters (i.e., J, A, R, Pmm, and 
Pmax), since these values are typically quite accurately known. 
However, it may not use values of the Cp-surface, the air 
density p, and the friction coefficient C/, since these values are 
inherently uncertain and can change over time. Furthermore, 
the controller should not rely on the wind speed Vw, since it 
may not be accurately measured. 

III. Controller Design 

In this section, we address the aforementioned problem 
by developing a nonlinear controller consisting of four sub- 
controllers. Figure [T] shows the architecture of the nonlinear 



controller, where each block represents a subcontroller. Note 
that the controller accepts Pd and Qd as reference inputs, 
uses i, ujr, P, and Q as feedback, and produces Vdr, Vqr, 
and (3 as control inputs to the wind turbine. Moreover, the 
different gray levels of the blocks in Figure [T] represent our 
intended time-scale separation in the closed-loop dynamics: 
the darker a block, the slower its dynamics. The subcontrollers 
will be described in Sections IIII-AljlirDl Note that Sections 
IIII-AI and IIII-BI up to the coordinate change are similar to our 
previous work ll25l . while the rest of the paper contains new, 
unpublished results. 

A. Rotor Voltages Subcontroller 

Observe that although the electrical dynamics ([T]i are non- 
linear, they possess a nice structure: the first and second 
rows of ^ are affine, consisting of linear terms and the 
constants Vds and Vqg, while the third and fourth are nonlinear, 
consisting of linear terms, the control variables Vdr and Vqr, 
and the nonlinearities —uirfgr and Urfdr induced by the 
electromechanical coupling. Since the nonlinearities enter the 
dynamics the same way the control variables Vdr and Vqr do, 
we may use feedback linearization [30J to cancel them and 
perform pole placement ISTl . i.e., let 



Vdr = i^r^gr — 
gr — ^r^dr 



- til, 



V, 



K2 'P + U2, 



(10) 

(11) 



where cUrPgr and —WrPdr are intended to cancel the nonlin- 
earities, —K'[ip and —K2f with Ki,K2 £ are for pole 
placement, and ui and U2 are new control variables to be 
designed in Section IIII-BI 

Substituting ( fTOb and ( fTTI ) into ([B, we get 



ip={A- BK)ip + [vd 



s Vqs Ul U2\ 



(12) 



where K — [Ki K2Y' the state feedback gain matrix. 
Since the electrical dynamics are physically allowed to be 
much faster than the mechanicals, we may choose K in (fT2l l 
to be such that A — BK is asymptotically stable with very 
fast eigenvalues. With K chosen as such and with relatively 
slow-varying ui and U2, the linear differential equation (fT2b 
may be approximated by a linear algebraic equation: 



ip ^ ~{A~ BK) ^ [vds Vqs Ul U2] 



(13) 



Consequently, the fifth-order state equations ([T]i and (|6]) may 
be approximated by the first-order state equation ^ along with 
algebraic relationships ( fTol i. (fTTT i. and ( fT3] l. This approximation 
will be made in all subsequent development (but not in 
simulation). 

Note that Q, ( fTOl l. and (fTTl i describe the Rotor Voltages 
Subcontroller block in Figure [U 

B. Electromagnetic Torque Subcontroller with Uncertainty Es- 
timation 

Having addressed the electrical dynamics, we now consider 
the mechanicals, where the goal is to construct a subcon- 
troller, which makes the rotor angular velocity track a 
desired, slow-varying reference Lo^d, despite not knowing the 
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Fig. 1. Model of the wind turbine and architecture of the nonlinear controller 



aerodynamic and mechanical parameters listed at the end of 
Section HI] 

To come up with such a subcontroller, we first introduce a 
coordinate change. As was shown in our previous work li25l . 
because of Q, dTji, and (ITSl l. the electromagnetic torque Tg 
may be expressed as a quadratic function of the new control 
variables ui and U2, i.e.. 



Te = [ui U2] 



11 
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where gi, q2, 93, 61, 62, and a depend on the electrical 
parameters and the state feedback gain matrix K. Moreover, 
as was shown in fZS], this quadratic function is always convex 
because its associated Hessian matrix [q\ ^3 ] is always positive 
definite. Since the mechanical dynamics in uj^, are driven 
by Tg, while Tg in ([T4l l is a quadratic function of ui and U2, 
the two new control variables ui and U2 collectively affect one 
state variable w^- This implies that there is a redundancy in ui 
and U2, which may be exploited elsewhere. Since the quadratic 
function is always convex, this redundancy may be exposed 
via the following coordinate change Ii25l . which transforms 
ui,M2 € M in a Cartesian coordinate system into r > and 
9 G [— tTjTt) in a polar coordinate system: 



2i 



6 = atan2(z2, zi) 



where 



2 



(15) 



(16) 



atan2() denotes the four-quadrant arctangent function, and M 
and D contain the eigenvectors and eigenvalues of [ %] on 
their columns and diagonal, respectively, i.e., A/-^ [W '^^]M — 
D. In the polar coordinates, it follows from ([T4li-(lT6]l that 

T,^r^+a', (17) 



where a' 



is always negative. From ^ and 
(ITtT i. we see that in the polar coordinates, is responsible 



for driving the mechanical dynamics in ujr and, hence, may 
be viewed as an equivalent electromagnetic torque, differed 
from Te only by a constant a'. On the other hand, the polar 
angle 9 has no impact on the mechanical dynamics and, 
thus, represents the redundancy that will be exploited later, 
in Section Hira 

Note that (ITSl l and dTSl l describe the Cartesian-to-Polar 
Coordinate Change block in Figure [T] 

Having introduced the coordinate change, we next show that 
the unknown aerodynamic and mechanical parameters, listed 
at the end of Section [III can be lumped into a scalar term, 
simplifying the problem. Combining (|6]l, dHJ, (|9]l, and (ITtT i. 



JiOr = 



-r^-a'-C 



U)r 



(18) 



Notice that the unknown parameters — namely, the Cp-surface, 
the air density p, the friction coefficient Cf, and the wind 
speed Via — all appear in (ITSl l. Moreover, these unknown 
parameters can be separated from the "control input" and 
lumped into a scalar function g{ujr, f3, Ko), defined as 



— a' — CfOJr 



(19) 



With g{Ldr, l3,Vi^) in ( lT9l ) representing the aggregated uncer- 
tainties, the first-order dynamics ( ITSl l are simpHfied to 



(20) 



To design a controller for r^, which allows the rotor angular 
velocity ujr to track a desired, slow-varying reference ujrd 
despite the unknown scalar function g{uJr, /?, V^), consider a 
first-order nonlinear system 

x^jif{x)+u), (21) 

where a; e M is the state, u e M is the input, and f{x) is a 
known function of x. Obviously, to drive x to some desired 
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value Xd E M, we may apply feedback linearization ||30l to 
cancel f{x) and insert linear dynamics, i.e., let 

u = -f{x) - a{x - Xd), (22) 

where a G M is the controller gain. Combining (ISTT i with 
yields the closed-loop dynamics 

a 

X = - Xd). 



(23) 



Thus, if a is positive, x in (l23l l asymptotically goes to Xd- 

Now suppose f{x) in (l2Tl i is unknown but a constant, 
denoted simply as / G M (we will relax the assumption that it 
is a constant shortly). With / being unknown, the controller 
(|22] | is no longer applicable. To overcome this limitation, 
we may first introduce a reduced-order estimator ||32i , which 
calculates an estimate / G K of /, and then replace f{x) in 
(l22T i by the estimate /: 



z = -j[u + f), 

f = z + hx, 

u = - f - a{x - Xd) 



(24) 

(25) 
(26) 



where z e K is the estimator state and h E R is the estimator 
gain. Defining the estimation error as / = f—f and combining 
(l2Tb with (|24]|-(|26T| yield closed-loop dynamics 

f = -j = -z-hx^--~f, (27) 



J" 

U - f - a{x - Xd)) 



-a{x-Xd)). (28) 



Hence, by letting both a and h be positive, both / and x in 
(|27] | and (l28T l asymptotically go to and Xd, respectively. 

Next, suppose both the state x and the desired value Xd 
must be positive, instead of being anywhere in M. With this 
restriction, the controller with uncertainy estimation (|24]|-(|26T| 
needs to be modified, because for some initial conditions, it is 
possible that x can become nonpositive. One way to modify 
the controller is to replace the linear term x — in (|26] l by 
a logarithmic one In — , resulting in 



Xd 



(29) 



With (|24] |. (l25T l. and ( |29] l, the closed-loop dynamics become 

(30) 



X = -{f ~a\n — ). 

J Xd 



(31) 



Note from (l3Tl i that for any / G M, there exists positive 
X, sufficiently small, such that x is positive. Therefore, for 
any initial condition (/(O), x(0)) with positive x(0), will 
remain positive, suggesting that the modification ( |29] | satisfies 
the restriction. 

Now suppose the input u must be nonpositive. With this 
additional restriction, (|29] | needs to be further modified. One 
way to do so is to force the right-hand side of ( |29] | to be 
nonpositive, leading to 



max{/ + a In — , 0}. 

Xd 



(32) 



Clearly, with ( [32] i, u is always nonpositive. 

Finally, suppose / is an unknown function of x, denoted 
as f{x). With this relaxation, we may associate the first-order 
nonlinear system (l2Tl i with the first-order dynamics (|20| | by 
viewing x as w^, Xd as cj^d, u as —r^, /(x) as ^(wr., /3, 14;) 
(treating /3 and 14j as constants), and f as g (i.e., 5 is an 
estimate of gitOr, /3, Vw)). Based on this association, ( |24| ). ( |25] ), 
and ( |32] | can be written as 

)i~r^+9), (33) 



J' 

g = z + huJr , 

r = metx{g + am ,0}. 



(34) 
(35) 



Having derived the controller with uncertainty estimation 
(l33]l-(l35Tl. we now analyze its behavior. To do so, some 
setup is needed: first, suppose uJrd, P, and are constants. 
Second, as was shown in ||25l , because of Assumptions (Al)- 
(A3) in Section [III there exists wf^-* E (0, 00) such that 
g(uj'^\ l3,Vw) = and g{ujr, l3,Vw) > for all ujr E 
iO,U^'^). Third, using ( fT9] l, and Assumptions (Al) and 
(A4), it is straightforward to show that there exist 7 E (—00, 0) 
and 7 E {0,oo) such that 7 < /?, 14,) < 7 for all 

Ur E (0, 00). Finally, with and (tj^-(l35Tl and with (w,-, g) 
as state variables (instead of (cj,., z)), the closed-loop dynamics 
can be expressed as 

\{g{ujr,f3,Viu) -ma,x{g + a\n—,0}), (36) 



g ^ z + hur 



J 



{g{uJr,l3,V^) - g). 



(37) 



The following theorem characterizes the stability properties 
of the closed-loop system ( |36] | and dJTl l: 

Theorem 1: Consider the closed-loop system ( 1361 ) and ( I37I ). 
Suppose UJrd, P, ond Vw are constants with < U)rd < ^r^\ 
where uj^r^ , along with 7 and 7, is as defined above. Let D = 

{{ujr,g)\0 < ujr < U!i'^\g E M.} C M^. If the controller gain 
a is positive and the estimator gain h is sufficiently large, i.e., 

(38) 



> 7 if 1> -|7, 



h > - 8(7+7) o'^^'^"^' 

f/ien.- (i) the system has a unique equilibrium point at 
{uJrdi gi^rd, l3,Vw)) in D; (ii) the set D is a positively 
invariant set, i.e., if (cUriO), g{0)) E D, then {ujr{t), g(t)) E D 
yt > 0; and (Hi) the equilibrium point {ujrd, g{^rd, 1^ iVw)) is 
locally asymptotically stable with a domain of attraction D. 
Proof See the Appendix. ■ 

Theorem [T] says that, by using the electromagnetic torque 
subcontroller with uncertainty estimation (l33Tl-(l35]l. if the 
gains a and h are positive and sufficiently large and if the 
desired reference iOrd does not exceed ujr , then the rotor 
angular velocity ujr asymptotically converges to w^d if uJrd, 
P, and Vw are constants and closely tracks ojrd if they are 
slow-varying. Notice that the gains a and h can be chosen 
independently of each other. Also, the condition "cOrd < Wr^^" 
is practically always satisfied, as uji^^ is extremely large 1251 . 

Note that ([33]l-(l35ll describe the Electromagnetic Torque 
Subcontroller with Uncertainty Estimation block in Figure [T] 
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Fig. 2. Relationships among the performance measure U, the to-be- 
determined variables 6, uj^d^ ™d (5, and the exogenous variables Vw, Pd< 
and Qjj. 



Equations (|40] i and ( 1411 1. which are represented in Figure |2] 
suggest that J7 is a function of the to-be-determined variables 
6, LUrd, and f3 as well as the exogenous variables P^, 
and Qd- Given that the smaller U is the better, these to-be- 
determined variables may be chosen to minimize U. However, 
such minimization is difficult to carry out because although Pd 
and are known, Vw is not. To make matter worse, since /i 
is known but /2 is not, the objective function is not entirely 
known. Somewhat fortunately, as was shown in Figure |2] 9 
affects U only through /i and not /2. Therefore, 9 may be 
chosen to minimize U for any given r^, uJrd, Pd, and Qd, i.e.. 



C. Polar Angle and Desired Rotor Angular Velocity Subcon- 
troller 

Up to this point in the paper, we have yet to specify how 9, 
UJrd, and (3 are determined. To do so, we first introduce a scalar 
performance measure and express this measure as a function 
of 9, UJrd, and /3. We then present a method for choosing these 
variables, which optimizes the measure. 

Recall that the ultimate goal is to make the active and 
reactive powers P and Q track some desired references Pd 
and Qd as closely as possible. Hence, it is useful to introduce 
a scalar performance measure, which characterizes how far P 
and Q are from Pd and Qd- One such measure, denoted as U, 
is given by 



9 = argmin^g[_^,,)/i(r^x,Wrd,-Pd,(9<i), 



(42) 



Wp Wpq 






Pd 


Wpq Wq 




Q- 


Qd_ 



(39) 



where Wp, Wq, and Wpq are design parameters satisfying Wp > 
and WpWq > Wpg, so that [ wf,' ] a positive definite 
matrix. With these design parameters, one may specify how 
the differences P — Pd and Q — Qd and their product {P — 
Pd){Q — Qd) are penalized. Moreover, with U in ( [39] l being 
a quadratic, positive definite function of P — Pd and Q — Qd, 
the smaller U is, the better the ultimate goal is achieved. 

Having defined the performance measure U via ( [39] l, we 
next establish the following statement: if the subcontrollers 
in Sections IIII-AI and IIII-BI are used with K chosen so that 
A — BK has very fast eigenvalues, a chosen to be positive, 
and h chosen to satisfy ( |38l l, and if 9, uJrd, P, Vw, Pd, and 
Qd are all constants, then after a short transient, U may be 
expressed as a known function /i of r^, 9, ujrd, Pd, and Qd, 
while r^, in turn, may be expressed as an unknown function 
/2 of Urd, P, and Vw, i.e.. 



U = hir\9,LUrd,Pd,Qd), 



(40) 
(41) 



as shown in Figure |2] To establish this statement, suppose the 
hypothesis is true. Then, after a short transient, it follows from 
(|39] | that t/ is a known function of P, Q, Pd, and Qd', from 
©, (II]!-©, ([Toll, and ([TT) that P and Q are known functions 
of (p, ojr. Ml, and U2', from (fT3] l that (p is a known function 
of ui and U2', from (flSl l and (fTSI l that ui and U2 are known 
functions of and 9; and from Theorem [T] that uj^ — uj^d- 
Thus, ( |40l i holds with /i being known. On the other hand, it 
follows from ( [35l l and Theorem [T] that j-^ = giuj^d, P^V^)- 
Hence, W\\ holds with /2 being unknown. 



which is implementable since r^, uj^d, Pd, and are all 
known. With 9 chosen as in ( |42] |. the minimization problem re- 
duces from a three-dimensional problem to a two-dimensional 
one, depending only on uj^d and /3. Since the objective function 
upon absorbing 9 is unknown and since Vw may change 
quickly, instead of minimizing U with respect to both ujrd 
and /3 — which may take a long time — we decide to sacrifice 
freedom for speed, minimizing U only with respect to ujrd 
and updating /3 in a relatively slower fashion, which will be 
described in Section UlI-DI 

The minimization of U with respect to uj^d is carried out 
based on a gradient-like approach as shown in Figure |3] To 
explain the rationale behind this approach, suppose /3, Vw, 
Pd, and Qd are constants. Then, according to (|40]|-(|42]|. U 
is an unknown function of ujrd- Because this function is not 
known, its gradient at any Wrd cannot be evaluated. To 
alleviate this issue, we evaluate U at two nearby w^d's, use the 
two evaluated [/'s to obtain an estimate of the gradient 
and move uj^d along the direction where U decreases, by an 
amount which depends on the gradient estimate. This idea is 
illustrated in Figure [3] and described precisely as follows: the 
desired rotor angular velocity uJrd{t) is set to an initial value 
ujrd{Q) at time t — Q and held constant until t — Ti, where 
Ti should be sufficiently large so that both the electrical and 
mechanical dynamics have a chance to reach steady-state, but 
not too large which causes the minimization to be too slow. 
From time t = Ti - Tq to t = Ti, the average of U{t), 
i.e., ifr Sxl-To U{t)dt, is recorded as the first value needed to 
obtain a gradient estimate. Similar to Ti, To should be large 
enough so that small fluctuations in U{t) (induced perhaps 
by a noisy Vw) are averaged out, but not too large which 
causes transient in the dynamics to be included. The variable 
uJrd{t) is then changed gradually in an S-shape manner from 
WrdCO) at time t = Ti to a nearby uJrd{^) + ^0Jrd{Ti) ait — 
Ti +T2, where Aujrd{Ti) is an initial stepsize, and T2 should 
be sufficiently large but not overly so, so that the transition 
in ujrd{t) is smooth and yet not too slow. The variable ujrd{t) 
is then held constant until t = 2Ti + T2, and the average 
of U{t) from t = 2Ti + T2 - To 10 t = 2Ti + T2, i.e., 
^ J'?Ti+T2-T U{t)dt, is recorded as the second value needed 
to obtain the gradient estimate. At time t — 2Ti + T2, the 
two recorded values are used to form the gradient estimate, 
which is in turn used to decide a new stepsize Aujrd{'^Ti+T2) 
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Fig. 3. A graphical illustration of the gradient-like approach. 



through 

AuJrd{2Ti+T2) 



ei sat 



_ 2T1+T2 

To J2T1+T2- 



U{t)dt-^J^l_^^U{t)dt\ 



(43) 



where ei > and £2 > are design parameters that define the 
new stepsize Aujrd{2Ti + T2), and sat() denotes the standard 
saturation function that limits Auirdi'^'Ti + T2) to ±ei. Upon 
deciding AuJrdi'^Ti + T2), 0Jrd{t) is again changed in an S- 
shape manner from Wr-d(O) + Awrd(Ti) at t = 2Ti + T2 to 
Wrd(O) + Aujrd{Ti) + Aa;rd(2Ti + T2) at t = 2Ti + 2T2, in 
a way similar to the time interval [Ti,Ti + T2]. The process 
then repeats with the second recorded value from the previous 
cycle [0, 2Ti + T2] becoming the first recorded value for the 
next cycle [Ti + T2,3Ti + 2T2], and so on. Therefore, with 
this gradient-like approach, ujrd is guaranteed to approach a 
local minimum when /3, Vw, Pd, and Qd are constants, and 
track a local minimum when they are slow-varying. 

Note that (|39]l, (|42]i, and (US describe the Polar Angle 
and Desired Rotor Angular Velocity Subcontroller block in 
Figure [T] 

D. Blade Pitch Angle Subcontroller 

As was mentioned, in order to speed up the minimization, 
we have decided to minimize U only with respect to LOrd, 
leaving the blade pitch angle /3 as the remaining undetermined 
variable. Given that an active power P that is larger than the 
rated value Prated of the turbine may cause damage, we decide 
to use /? to prevent P from exceeding Prated, thereby protecting 
the turbine. Specifically, we let /3 be updated according to 



^e3(^rated ^ P) 



if /3 = /3, 
if /3 = /3, 
otherwise, 



and P < Prated, 

max 

and P > Prated, 



(44) 



where £3 > is a design parameter that dictates the rate at 
which (3 changes. Note that with ( |44] |. (3 is guaranteed to lie 
between /3min and /3max- Moreover, when P is above (below) 
P increases (decreases) if possible, in order to ti-y to 
capture less (more) wind power, which leads to a smaller 
(larger) P. 

Note that (l44l l describes the Blade Pitch Angle Subcontroller 
block in Figure [T] 



Remark 1: The blade pitch angle subcontroller may be 
designed based on other considerations. For example, if the 
forecast of, say, the hourly-average wind speed Vw is avail- 
able, for blade protection (3 may be chosen as /? = F{Vw) 
for some non-decreasing function F : (0, 00) — )■ [/3niin, Anax] 



IV. Simulation Results 

To demonstrate the capability and effectiveness of the pro- 
posed controller, simulation has been carried out in MATLAB. 
To describe the simulation settings and results, both the per- 
unit and physical unit systems will be used interchangeably. 

The simulation settings are as follows: we consider a 
1.5 MW, 575 V, 60 Hz wind turbine that is essentially adopted 
from the Distributed Resources Library in MATLAB/Simulink 
R2007a. The values of the wind turbine parameters are: 
= Ipu, Rs = 0.00706 pu, Rr = 0.005 pu, ^ 3.071 pu, 
Lr — 3.056 pu, Lm = 2.9 pu, Vds — 1 pu, Vqs — Opu, 
J = 10.08 pu, A = 4656.6 m^, R = 38.5 m, P^-.^ = Odeg, 
/^max — 30deg, and C/ = 0.01 pu. The Cp-surface adopted 
by MATLAB, which is taken from 133], is Cp(A,^) = 

ci - C3/3 - C4) e-^ + ceA, where J- = ^^^^ - 
ci = 0.5176, C2 = 116, C3 = 0.4, a = 5, C5 = 21, 
and C6 = 0.0068. The mechanical power captured by the 
wind turbine is P„i(pu) = p^"" Cp(pu) Vu,(pu)'^, where 
P„j(pu) — -p^, Pnom = 1.5 MW is the nominal mechan- 
ical power, Pwind_ba,se ~ 0.73 pu is the maximum power at 
the base wind speed, Peiec_base — 1-5 x 10^/0.9 VA is the 
base power of the electrical generator, Cp(pu) 



c. 



v: 



p_nom = 0.48 is the peak of the Cp-surface, T4,(pu) — y ^ , 
and yu)_base = 12m/s is the base wind speed. Note that' the 
maximum mechanical power, captured at the base wind speed. 



r(pu) 



is 0.657 pu. The tip speed ratio is A(pu) 

'*nnm 



where 



V,„(pu)' 

A(pu) = Anom = 8.1 is the A that yields the peak of 

the Cp-surface, Wr_base = 1.2 pu is the base rotational speed, 
Wr(pu) = , and ojr nom = 2. 1039 rad/scc is the nominal 
rotor angular velocity. For more details on these parameters 
and values, see the MATLAB documentation. 

As for the proposed controller, we choose its parame- 
ters as follows: for the Rotor Voltages Subcontroller, we 
let the desired closed-loop eigenvalues of the electrical 
dynamics be at -10, -15, and -20 ± 5j. Using MAT- 
LAB 's place function, the state feedback gain matrix 
K ~ [Ki K2Y' that yields these eigenvalues is found 
to be Ki = [12277 -4493.8 32.7 -6.3]^ and K2 = 
[-1615.4 12117 -0.4 32.2]^. Moreover, we let a = 5 and 
h — \Q for the Electromagnetic Torque Subcontroller with 
Uncertainty Estimation; Wp = 10, Wq — 1, Wpq = 0, 
ei = 0.025, 62 = 2, To = Is, Ti = 4s, and T2 = 6s 
for the Polar Angle and Desired Rotor Angular Velocity 
Subcontroller; and 63 = 2.7 and Prated = 1 pu for the Blade 
Pitch Angle Subcontroller. 

The simulation results are as follows: we consider a scenario 
where the wind speed V^ is derived from an actual wind 
profile from a wind farm located in northwest Oklahoma, 
the desired active power Pd experiences large step changes 
between 0.3 pu and 1 pu, and the desired reactive power Qd 
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Fig. 4. Effective operation in both the MPT and PR modes and seamless switching between them under an actual wind profile from a wind farm located in 
northwest Oklahoma. 



is such that the desired power factor PF^ is fixed at 0.995. 
As will be explained below, these values force the turbine to 
operate in both the MPT and PR modes, along with switching 
between them, under a realistic wind profile. Figure |4] shows 
the simulation results for this scenario in both the per-unit 
and physical unit systems, wherever applicable. Subplot 1 
shows the wind profile Vw Subplot 2 shows the value of 
Cp, while subplot 3 shows the desired and actual active 
powers Pd and P. Note that, for the first 1200 seconds during 
which Pd is unachievable at 1 pu, the turbine operates in the 
MPT mode and maximizes P, as indicated by the value of 
Cp approaching its maximum of 0.48 after a short transient 
(the turbine is initially at rest). At time 1200s when Pd 
drops sharply from 1 pti to an achievable value of 0.3 pu, the 
turbine quickly reduces the value of Cp, accurately regulates 
P around Pd, and effectively rejects the "disturbance" V^, 
thereby smoothly switches from the MPT mode to the PR 
mode. At time 2400s when Pd goes from 0.3 pti back to 1 pu, 
the MPT mode resumes. Because Vw is strong enough at that 
time, P approaches Pd- Moreover, the moment P exceeds Pd 
(which is equal to Prated), the blade pitch angle (3 increases as 
shown in subplot 8 in order to clip the power and protect the 



turbine. At time 3275s when Vw becomes weaker, j3 returns to 
/3min = deg, thereby allowing the value of Cp to return to its 
maximum of 0.48 and P to be maximized. Subplot 4 shows the 
desired and actual power factors PF^ and PF, while subplots 5 
and 6 show the desired and actual rotor angular velocities 
ujrd and Ur in normal and zoomed-in views. As can be seen 
from these subplots, throughout the simulation, both PF and 
ujr are maintained near PF^ and Wrd, respectively, affected 
only slightly by the random wind fluctuations. Moreover, the 
small S-shape variations in Urd in subplot 6 resemble those in 
Figure |3] Finally, subplots 7 and 8 show the control variables, 
i.e., the rotor voltages Vdr and Vqr as well as the blade pitch 
angle f3. 

The above simulation results suggest that the proposed con- 
troller not only is capable of operating effectively in both the 
MPT and PR modes, it is also capable of switching smoothly 
between them — all while not knowing the Cp-surface, air 
density, friction coefficient, and wind speed. 

V. Concluding Remarks 

In this paper, we have designed a controller for a variable- 
speed wind turbine with a DFIG. The controller, consisting 
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of four subcontrollers, has been developed based on a fifth- 
order, electromechanically-coupled, nonhnear model of the 
wind turbine by integrating several linear and nonlinear control 
strategies and exploiting time-scale separation in the dynamics. 
We have shown that the controller is able to make the wind 
turbine operate in both the MPT and PR modes and switch 
smoothly between them, in addition to maintaining a desired 
power factor. Furthermore, the controller does not require 
knowledge of the Cp-surface, air density, friction coefficient, 
and wind speed. Simulation has been carried out using a 
realistic wind profile, and the results demonstrate the capability 
and effectiveness of the controller Future work includes de- 
signing a comprehensive wind farm controller for a multitude 
of turbines, which builds upon results documented in this 
paper. 

Appendix 
Proof of Theorem[T] 

First, we show (i). Setting Ur and g in ( [36l l and (l37l l 
to zero yields g{ijjr, l3,Vw) = max{g + a\n -^,0} and 
g = g{uJr, P, Vui)- When g + a\n > 0, we have ujr — ^rd 
and g = K;)- Thus, (ujrd, g{uJrd, l3,Vw)) is an 

equilibrium point, which is in D, since < iVrd < uji^\ 
On the other hand, when q + a\n < 0, we have e 
and g — 0, where i7 = {oj e (0,oo) : g(a;,/3, — 0} and 
ul'^'' — min n. Since g + aln < and g = 0, we have 

ujr < i^rd- Since ujr € i^r^'' = mln 17, and ujrd < i^r^\ we 
have ojr > oJrd- Hence, there is a contradiction, implying that 
when g + aln-^ < 0, there is no equilibrium point in D. 
This proves (i). 

Next, we show (ii). To do so, it is useful to think of D 
as a vertical strip in the two-dimensional state space {ujr,g). 
Notice that on the right boundary of the strip where ujr = 



Differentiating V and using ( |36] | and i37\ . we get 



,(1) 



(jjr ', because of (l36l l and because g{uj'!.^\ l3,Vw) = and 
max{g + aln^fy,0} > 0, we have Ur < 0. Thus, the state 
{uJr,g) cannot escape D through the right boundary. Next, 
note that for each fixed g £ M., there exists w* > such that 
for all Ur e (0,a;*), .g + aln ^ < 0. This, along with ^ 

and the fact that g{uj, /3, Vw) > for all u e (0, Wr^^), implies 
that near the left boundary of the strip where ujr is arbitrarily 
small but positive, we have ujr > 0. Hence, the state {u!r,g) 
cannot escape D through the left boundary. This proves (ii). 

Finally, we show (iii). Consider a Lyapunov function candi- 
date V : D ^ M., defined as V{u!r, g) = ac{uJr In — cj^ + 
i^rd) + ^{g{i^r, /?, Vw) ^g)'^, where c > is to be determined. 
Note that V is continuously differentiable over D. Moreover, 
V is positive definite over D with respect to the equilibrium 
point {ujrd,g{uJrd,P,Vn,)), since V{uJrd, g{uJrd, ^,Vn,)) = 
and V{uJr,g) > for all {ujr,g) ^ {ujrd, g{uJrd, P^V^n j) due 
to the property ojr In ^ — ujr + ^rd > for all iVj. ^ uj^d- 
Furthermore, V is unbounded toward the top, bottom, and 
left boundary of the vertical strip D, but not so toward 
the right boundary of D. This is because for each fixed 
uJr G (OjWr^''], lim|g|^oo g) — oo, and for each fixed 

5 e M, lim^^^^o ^(<^ri ff) = oo and V{Lijl^\g) < oo. Note 
that although V is not unbounded toward the right boundary 
of D, the state {ujr,g) cannot cross this boundary due to (ii). 



JV= 



^{g-g) 



max{g + aln^,Oy 

Hg - g) 



where, for convenience, the function arguments are omitted. 
Note that because of (ii) and the above properties of V, to 
show (iii), it suffices to show that V is negative definite over 
D with respect to the equilibrium point {ujrd, gi^rd, P ,Vw))- 
To this end, let D be partitioned into two disjoint sets 
Di = {{ujr.g) e Z? : .g + a In ^ > 0} and Da = {K, g) & 
D : g + a\n < 0}. Note that the equilibrium point 
{ujrd,g{i^rd,f3,Vj)) is in Di. 

Suppose {uJr,g) E Di. Then, V takes a quadratic form: 



JV 



"In 


T ' 


<^rd 




9-g. 
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a c 
2 va^r 



c) 



a ( dg 
2 ^duir- 



c) 



A. 



In^ 
.9 - .9 



Note that if {iOr,g) = {^rd, gi^^rd, P,Vw)), = 0. Also, the 
leading principal minors of the above symmetric matrix are 
a^c and a'^c{h — ^(^j- + c)^). Thus, if h and c satisfy 



C)2 >0, 



Vu;e(0,oo), (45) 



then this symmetric matrix is positive definite, so that V < 
for any {ujr,g) {^rd, g{<^rd, f3,Vw))- Therefore, if h and c 
satisfy ( |45] ). V is negative definite over Di with respect to 

{uJrd,giuJrd,l3,Vw))- 

Next, suppose (w,., g) e D2- Then, V is bounded from 
above by a quadratic form: 



JV = -hg^ 
< -hg^ 





T ' 


g 




_g_ 





i2h-^)ggH^ 

OUJr OUJr 

i2h-^)ggH^ 

OUJr OUJr 



h 

2 V aw, ^ ^ 



h)g' 
h)g' 



2h) 



l( dg 

2^aw, 



UJr 

acg m 

i^rd 



egg 

2h) 



OUJr 









g 




_g_ 



Note that the leading principal minors of the above symmetric 
matrix are h and c{h — + c)^)- Thus, if h and c 

satisfy ( |45l l, then this symmetric matrix is positive definite. 
Since (uJr^g) G D2 and ujrd < uji^\ if g = 0, then g > 0. 
Thus, g and g cannot be zero simultaneously. Hence, V < 0. 
Therefore, if h and c satisfy i45[ . V is negative over D2- 

As it follows from the above, if h and c satisfy (l45T l. V is 
negative definite over D with respect to the equilibrium point 
{ujrd, gi^rd, f3,Vw)), SO that (iii) holds. 

It remains to show that if h satisfies dSSl l. then there 
exists c > such that (l45l l holds. Suppose h satisfies dSST l. 

Let F{2,j) = 7 if 7 > -^7 and ^(7,7) -^^7^ 
otherwise. Then, h > ^(7,7). Let f{x,j,j) — max{(7 + 



^(7 



■}, where a; > 0. Then, it can be shown that 



F{"f,^) — minj;>o /(a:, 7, 7) by considering the following 



and 



three cases separately: 7 > —7, —7 > 7 > —57, 
— hi > 7- Because h > F{j,^), there exists c > 0, given by 



c = argmina:>o /(x, 7, 7), such that h > /(c, 7,7). Because 
7 ^ '^di'^' ^) — 7 f^i" all w S (0, 00) and by definition 
of /(a;, 7, 7), we have ^(£5(0;, ^, K,) + c)^ < /(c,7,7) for 
all uj £ {0,00). Since h > /(c, 7,7), ( l45T l holds, as desired. 
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